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Typical Pulsar Analysis 

► Download data and screen events 



-Event file 
-Spacecraft file 

-Pulsar database file (Fermi_PulsarDB_v001 .fits) 


► Analyze time series data 

-Search for pulsations 

-Determine or refine pulsar ephemeris 

-Confirm periodicity published in the literature 



► Assign pulse phase or orbital phase to each photon 

-Folded light curve (pulse shape or pulse profile) 

-Phase-resolved spectrum 
-Phase-resolved image 
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Tips for Data Screening 


►For time series analysis 

- For maximum sensitivity on pulse detection and frequency 
determination, screen data to maximize pulsed component and 
minimize constant (“background”) component. 

- Tight spacial selection (typically size of a couple of PSFs) 


- Wider energy range 

- Broader choice of event types 


►For analysis using pulse/orbital phase 

- Use standard event selection for spectral or image analysis. 
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Time Series Analysis 

► Periodicity test 

- Chi-squared test (Leahy et at. 1983, ApJ 266, 160; Vaughan et al. 1994, 
ApJ 435, 362) 

- Z n 2 test (Buccheri et al. 1983, A&A 128, 245) 

- Rayleigh test (equivalent to Z n 2 test for n = 1) 

- H test (De Jageret al. 1989, A&A 221, 180) 

• Science Tools: gtpsearch, gtptest 

• Other tools: ef search (XRONOS) 

► Period search 

- Discrete Fourier transforms (Ransom 2002, AJ 124, 1788) 

• Science Tools: gtpspec 

• Other tools: powspec (XRONOS) 
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Pulse / Orbital Phase Assignment 


>Pulse Phase 

- Computes a pulse phase for each photon. 

- Appends the phase value to the event entry for later use. 

• Science Tools: gtpphase 

• Other tools: TEMP02 with Fermi plugin 


► Orbital Phase 

- Computes an orbital phase for each photon. 

- Appends the phase value to the event entry for later use. 

• Science Tools: gtophase 
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Utility Tools 

»Ephemeris computer gtephem 





- Reads pulsar ephemerides database and computes pulsar’s spin 
ephemeris (such as pulse frequency) at a given moment in time. 

- Also serves as a sanity checker for first-time users of our pulsar 
ephemerides database. 


► Pulsar ephemerides manipulation tool gtpulsardb 

- Sub-selects pulsar ephemerides from a master database. 

- Also used to create, modify, merge pulsar ephemerides database. 

► Photon arrival time correction gtbary 

- Computes barycentric or geocentric time for each event and modifies 
event time in event file, for use with non-Fermi analysis tools such as 
XRONOS tools. 

CAUTION: event files processed by gtbary must not be used in any 
analysis with Science Tools. 
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Analyzing Vela data 


prompt)- gtpspec 

This is gtpspec version Sc ienceToo is -v9rl5p2-fssc -20090701 
Event data file name[] Vela_lG0_3QQQQ0_evtQ2*f its 
Spacecraft data file name[] L090923110451E0D2F37E17_SC00*f its 

Name of input pulsar ephemerides database file (for binary demodulation only)[] Fermi_PulsarDB_vQ01*fits 

Pulsar name (for binary demodulation only)[RNV] PSR J0835-4510 

Output FITS file name (NONE for no FITS output) [] gtpspec_out + f its 

Uidth of time bins, in seconds (0 * i ) [l*e-2] 0*01 

Number of time bins to be transformed at once [1000000] 10000000 

How will the time origin of the periodicity test be specified? (STRRT | STOP | MIDDLE | USER) [MIDDLE] MIDDLE 
Right Rscension to be used for barycenter corrections (degrees) [0 * ] 128*836048 
Declination to be used for barycenter corrections (degrees) [0 + ] -45*176425 
How will spin ephemeris be specified? (FREQ | PER) [FREQ] FREQ 

Ratio of frequency first time derivative to frequency at the time origin (Hz)[0*] 0*0 
Ratio of frequency second time derivative to frequency at the time origin (Hz/s)[0*] 0*0 
Search Type: Fourier Rnalysis 
Fourier Resolution: le-05 Hz 
Sampling Frequency: le-05 Hz 

Data Binning: 27 segments with 10000000 time bins in each segment 
Probability Distribution: Chi Squared with 54 degrees of freedom 
Search Range (Hz): [0*01, 50] 

Number of Trial Frequencies: 4999000 

Number of Independent Trials: 4999000 

Maximum Statistic: 518*813411319925 at 11*19001 Hz 

Chance Probability Range: (1 *751880806 14015e -70, l*75205599227056e-70) 

Warning in <TClass: : "DC lass): no dictionary for class st_graph::IFrame is available 
prompt) 
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Analyzing Vela data (cont.) 


Fourier Analysis: Power Spectrum 


500 

400 

V- 

oj 

3 

£300 

1 — 

s 

£200 

100 

0 



Pulse Frequency (Hi) 
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Analyzing Vela data (cont.) 


prompt> gtpsearch 

This is gtpsearch version ScienceTools-v9rl5p2-fssc-20090701 
Event data file name[] Vela_100_300000_evt02.f its 
Spacecraft data file name[] L090923 1 1045 1E0D2F37E17_SC00. fits 
Pulsar ephemerides database file name[] Fermi_PulsarDB_v001.f its 
Pulsar name [ANY] PSR J0835-4510 

Output FITS file name (NONE for no FITS output) [] gtpsearch_out.fits 

Type of statistical test to perform (CHI2 - Chi squared, RAYLEIGH - Rayleigh test, Z2N - Z2n test, H - H test) (CHI2 | RAYLEIGH |Z2N |H) [CHI2] Chi2 
Number of phase bins for Chi2 test [10] 10 

Size of steps for trials, in units of the Fourier resolution (0.:) [0.5] 0.5 
Number of trials [100] 600 

How will the time origin of the periodicity test be specified? (START | STOP | MIDDLE | USER) [MIDDLE] USER 
Time origin of the periodicity test[0.] 55028.5 

Time format for the user -supplied time origin (FILE | MJD | ISO | FERMI | GLAST) [FILE] MJD 
Time system for the user -supplied time origin (FILE j TAI j TDB j TT | UTC) [FILE] TDB 
How will spin ephemeris be specified? (DB | FREQ | PER) [DB] FREQ 
Epoch for the spin ephemeris[0.] 54663.0 

Time format for spin ephemeris epoch (FILE | MJD | ISO | FERMI | GLAST) [FILE] MJD 

Time system for spin ephemeris epoch (FILE j TAI j TDB j TT | UTC) [FILE] TDB 

Right Ascension to be used for barycenter corrections (degrees) [0.] 128.836048 

Declination to be used for barycenter corrections (degrees) [0.] -45.176425 

Pulse frequency at the epoch of the spin ephemeris (Hz) (0.:) [1.] 11.19051540397055 

First time derivative of the pulse frequency at the epoch of the spin ephemeris (Hz/s)[0.] - 1.559072535 133380e-ll 

Second time derivative of the pulse frequency at the epoch of the spin ephemeris (Hz/s/s)[0.] 0.0 

Search Type: Folding Analysis 

Fourier Resolution: 3.73754e-07 Hz 

Sampling Frequency: 1.86877e-07 Hz 

Type of test: Chi -squared Test, 10 phase bins 

Probability distribution: Chi -squared, 9 degrees of freedom 

Search Range (Hz): [11.1899669982264, 11.1900789375696] 

Number of Trial Frequencies: 600 
Number of Independent Trials: 300 

Maximum Statistic: 3617.22319536 at 11.1900234350906 Hz 


Chance Probability Range: (0, 6. 11271 140709 163e -97) 

Uarning in CTClass: :TClass>: no dictionary for class st_graph: :IFrame is available 
prompt> 
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Analyzing Vela data (cont.) 
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Analyzing Vela data (cont.) 



pro mpt > gtpphase 

This is gtpphase version Sc ienceToo Is -v9rl5p2-fssc -20090701 
Event data file name[] Vela_lQ0_300000_evt02 + f its 
Spacecraft data file name[] L090923110451E0D2F37E17_SC00 + f its 
Pulsar ephemerides database file nadie[] Fermi_PulsarDB_v001*fits 
Pulsar name[ANY] PSR J0835-4510 

How will spin ephemeris be specified? (DB | FREQ | PER) [DB] FREQ 
Epoch for the spin ephemeris [0 + ] 55028*5 

Time format for spin ephemeris epoch (FILE | MJD | ISO | FERMI | GLRST) [FILE] MJD 
Time system for spin ephemeris epoch (FILE | TRI |7DB |TT|UTC) [FILE] TDB 
Right Rscension to be used for bar y center corrections (degrees) [0*] 128*836048 
Declination to be used for bar y center corrections (degrees) [0 * ] -45*176425 
Base value of phase at this epoch [0*] 0*5 

Pulse frequency at the epoch of the spin ephemeris (Hz) (0*:) [1*] 11*1900234350 906 

First time derivative of the pulse frequency at the epoch of the spin ephemeris (Hz/s)[0 t ] -l t 559072535133380e-ll 

Second time derivative of the pulse frequency at the epoch of the spin ephemeris (Hz/s/s) [0 + ] 0 + 0 

prompt> 
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Analyzing Vela data (cont.) 



prompt> fhisto louva 1=0*0 highvai=l + 0 
Name of FITS file and [ext#][] Vela_100_300000_evt02 + f its 
Name of output FITS file [] velapulsar_pulse_shape*f its 
Name of column to generate histogram [] PULSE_PHflSE 
Size of bins[INDEF] 0*05 

input file (no b inspec ) : file: // Ve la_ 100_300000_e vt02 * f it s [ 1 ] 

input file : f ile ://Vela_100_300000_e vt02 *f its [1] [bin PULSE_PHRSE=0 : 1:0 *05] 

The histogram extension: ldhisto 

prompt> 


prompt> fplot 

Name of FITS file and [ext#][] velapulsar_pulse_shape*f its 
Name of X fix is Parameter [error] [] X 

Name of V fixis Parameter [error] up to S allowed [] V [Error] 
Lists of rows[-] - 

Device: /X Window, /XTerm, ,/TK, /PS, etc [/XU] /xserv 
finy legal PLT command [] marker on 
PLT> exit 
prompt> 
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Analyzing Vela data (cont.) 


Plot of file velapulsar_puls«_shope.fits 



29 -Sep- 2009 16:26 
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Analyzing Vela data (cont.) 


prompt> gtptest 

This is gtptest version Sc ienceToo Is -v9rl5p2-fssc -20090701 

Event data file name[] Vela_100_300000_evt02.f its 

Output FITS file name (NONE for no FITS output) [] gtptest _out .f its 

Number of phase bins for Chi2 test [10] 10 

Number of harmonics for Z2n test [10] 10 

Maximum number of harmonics for H test [10] 10 

Type of test: Chi-squared Test, 10 phase bins 

Probability distribution: Chi -squared, 9 degrees of freedom 

Test Statistic: 3617.22319536 

Chance Probability Range: (0, 2.03757046903054e-99) 

Type of test: Rayleigh Test 

Probability distribution: Chi -squared, 2 degrees of freedom 


Test Statistic: 1700.6017260263 

Chance Probability Range: (0, 2.00883927982452e-99) 

Type of test: Z2n Test, 10 harmonics 

Probability distribution: Chi -squared, 20 degrees of freedom 

Test Statistic: 4797.77337691268 

Chance Probability Range: (0, 2.07853386442652e-99) 

Type of test: H Test, 10 maximum harmonics 
Probability distribution: H Test -specific 
Test Statistic: 4762.16995611264 



Chance Probability Range: (0, 4e-08) 

Warning in <TC lass : : "DC lass > : no dictionary for class st_graph: :IFrame is available 
prompt> 
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Analyzing Vela data (cont.) 

Folded! Light Curve 
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